Version date: December 16, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



THE VLA SURVEY OF THE CHANDRA DEEP FIELD SOUTH. 
IV. SOURCE POPULATION 

P. Padovani, V. Mainieri 

European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei Miinchen, Germany 

P. TOZZI 

INAF, Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 1-34131, Trieste, Italy 

K. I. Kellermann, E. B. Fomalont 

National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903-2475, USA 

N. Miller 

Department of Physics and Astronomy, Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218, USA 

P. Rosati, P. Shaver 

European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei Miinchen, Germany 

Version date: December 16, 2008 

ABSTRACT 

We present a detailed analysis of 256 radio sources from our deep (flux density limit of 42 /jjy 
at the field centre at 1.4 GHz) Chandra Deep Field South 1.4 and 5 GHz VLA survey. The radio 
population is studied by using a wealth of multi-wavelength information in the radio, optical, and 
X-ray bands. The availability of redshifts for ~ 80% of the sources in our complete sample allows us 
to derive reliable luminosity estimates for the majority of the objects. X-ray data, including upper 
limits, for all our sources turn out to be a key factor in establishing the nature of faint radio sources. 
Due to the faint optical levels probed by this study, we have uncovered a population of distant Active 
Galactic Nuclei (AGN) systematically missing from many previous studies of sub-millijansky radio 
source identifications. We find that, while the well-known flattening of the radio number counts below 
1 mjy is mostly due to star forming galaxies, these sources and AGN make up an approximately 
equal fraction of the sub-millijansky sky, contrary to some previous results. The AGN include radio 
galaxies, mostly of the low-power, Fanaroff-Riley I type, and a significant radio-quiet component, 
which amounts to approximately one fifth of the total sample. The ratio of radio to optical luminosity 
depends more on radio luminosity, rather than being due to optical absorption. 

Subject headings: radio continuum: galaxies — X-rays: galaxies — galaxies: active — galaxies: star- 
burst 



1. INTRODUCTION 

The extragalactic radio source population ranges from 
normal galaxies with luminosities ~ 10 20 W Hz -1 to 
galaxies whose radio emission is as much as 10 4 — 10 7 
times greater owing to regions of massive star formation 
or to an active galactic nucleus (AGN). The population 
of radio sources in the sky with flux densities > 1 mjy 
is dominated by AGN driven emission generated from 
the gravitational potential associated with a supermas- 
sive black-hole in the nucleus. For these sources, the 
observed radio emission includes the classical extended 
jet and double lobe radio sources as well as compact ra- 
dio components more directly associated with the energy 
generation and collimation near the central engine. Be- 
low 1 mjy there is an increasing contribution to the radio 
source population from synchrotron emission resulting 
from relativistic plasma ejected from supernovae associ- 
ated with massive star formation in galaxies or groups 
of galaxies, often associated with mergers or interac- 
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tions (e.g.. IWindhorst et all 119951: iRichards et "all 119981 : 
iFomalont et al.ll2002ft .' However, the mix of star forming 
galaxies (SFG) and AGN and their dependence on epoch 
is not well dete rmined. 

For example. iGruppioni et al.1 (|2003l ). by using 15-^m 
Infrared Space Observatory (ISO) and 1.4 GHz data 
have suggested that starburst galaxies start to become 
a significant fraction of the radio source population for 
S T < 0.5-0.8 mjy, e xceeding ~ 60% only at S r < 50 //Jy. 
On the other hand, iHuynh et al.l ((2005) , by fitting ra- 
dio number counts of the Hubble Deep Field South and 
other surveys with evolutionary starburst models, have 
estimated that the number of SFG exceeds 50% of the 
total already at ~ 0.25 mjy. 

The most commonly accepted paradigm has been that 
the sub-mJy population is largely made up of SFG. A 
(small but growing) number of researchers, however (e.g., 
Gruppioni et al.lll99 9: Cilie gi et al.ll2003t ISevmour et al.l 
20081 : ISmolcic et al.l 120081 ) . have pointed out that radio 
sources (radio galaxies and radio-loud AGN) could also 
make a significant contribution to the sub-mJy popula- 



2 



Padovani et al. 



tion It has been further suggested (|Jarvis fc Rawlingsl 
12004 iSimpson et al.ll20M ) that radio-quiet AGN (see be- 
low) might become relevant at sub-mJy levels. 

The Chandra Deep Field South (CDFS) area is part of 
the Great Observatories Origins Deep Survey (GOODS) 
and as such is one of the most intensely studied region 
of the sky. High sens itivity X-ray observa t ions are avail- 
able f rom Chandra ( | Giacconi et all 120021 : iLehmer et al.l 
[20051: ILuo et all 12008ft . The GOODS and Galaxy Evo- 
lution from Morphology and Spectral energy distribu- 
tions (GEMS) multiband imaging programs using the 
HST Advanced Camera for Surveys (ACS) give sen- 
sitive high resolution optical images (|Giavalisco et al.l 
[200l:lRixet"aII[200l . Ground based imaging and spec- 
troscopy are available fr om the ESQ 2.2m a nd 8m tele- 
scopes (see references in lMainieri et al. 2008, Paper II). 

We observed the CDFS with the NRAO Very Large 
Array (VLA) for 50 hours at 1.4 GHz mostly in the BnA 
configurations in October 1999 and February 2001, and 
for 32 hours at 5 GHz mostly in the C and CnB con- 
figuratio ns between June and O ctober 2001 (details are 
given m lKellermann~ al. 2008, Paper I). The effective 
angular resolution was 3.5" and the minimum root-mean- 
square (rms) noise was as low as 8.5 /iJy per beam at 1.4 
GHz and 7 ^tJy per beam at 5 GHz. A total of 266 radio 
sources were catalogued at 1.4 GHz, 198 of which are in a 
complete sample with signal-to-noise ratio (SNR) greater 
than 5 and located within 15' of the field center. The cor- 
responding flux density ranges from 42 /zJy at the field 
center to 125 ^tJy near the field edge. These deep radio 
observations complement the larger area, but less sensi- 
tiv e, lower resolution o bservati ons of the CDFS disc ussed 
bv lAfonso et~ail (|2OO60 and bv lNorris et~ai] (|2006fh 

In addition to our deep radio data, our set of ancil- 
lary data is quite unique and allows us to shed new light 
on the nature of the sub-mJy radio source population. 
These data include reliable optical/near-IR identifica- 
tions for 94% of the radio sources, optical morphological 
classification for 61% of the sample, redshift information 
for 73% of the objects, and X-ray detections for 33% of 
the sources, with upper limits for all but three of the 
others. This is the largest and most complete sample of 
/iJy sources in terms of redshift information. 

The purpose of this paper is to study the composi- 
tion of the radio-detected VLA-CDFS sample in terms 
of source properties by using all available information. 
More details on the VLA-CDFS can be found in Paper 
I, which addresses the radio data, in Paper II, which 
discusses the optical and near IR counterparts to the ob- 
served radio sources, and in iTozzi et ail (|2008fl (Paper 
III), which studies the X-ray properties. 

The sample radio-optical data are described in § 2, 
while § 3 deals with the X-ray data. § 4 analyzes the 
overall properties of our sample and derives the radio 
number counts for various classes together with some 
model predictions, while § 5 discusses our results. Fi- 
nally, § 6 summarizes our conclusions. Throughout this 
paper spectral indices are defined by S v oc v~ a and the 
values H = 70 km s _1 Mpc -1 , = 0.3, and £1a = 0.7 
have been used. Th is paper supersedes the preliminary 
results presented bv lPadovani et all (|2007bf ). 



2. THE SAMPLE RADIO-OPTICAL PROPERTIES 



2.1. Expected radio properties 

Before we turn to the radio sources, we first define 
the typical ranges in radio power and radio-to-optical 
flux density ratio for the type of objects which are most 
likely to be part of our survey: galaxies (normal, radio, 
and star forming), and active galactic nuclei (AGN). 

Normal galaxies, with radio powers < 4 x 10 W 
Hz -1 , are thought to have their radio emission dom- 
inated by synchr otron radiation from interstellar rela - 
tivistic electrons ([Philips et al.ll 19861 : [Sadler et al.lll989h . 
Radio galaxies are considered to have gigahertz ra- 
dio power > 10 22 W Hz" 1 (e.g., ISadler et al.l fl989l: 
iLedlow fe Owenlll996| ). associated with relativistic jets, 
kiloparsec scale morphology, which extends well beyond 
the region of optical luminosity, and, at least at low 
redshift, mostly elliptical hosts. This radio power is at 
the faint end of the luminosity f unction of radio galax- 
ies (e.g.. lUrrv fc Padovanil [i~995l ) and, therefore, repre- 
sents a 'natural' threshold for " radio-loudness" in galax- 
ies, although recent papers have suggested the pres- 
ence of non-thermal, jet emissio n in early type galaxies 
as fai nt as 10 20 W Hz -1 : e.g., iBalmaverde fc Capettil 
( 2006). iFanaroff fc Rilevl (|1974f) recognized that radio 
galaxies separate into two distinct luminosity classes, 
each with its own characteristic radio morphology. High- 
luminosity Fanaroff-Riley (FR) lis have radio lobes with 
prominent hot spots and bright outer edges, while in 
low-luminosity FR Is the radio emission is more diffuse. 
The luminosity distinction is fairly sharp at 178 MHz, 
with FR Is and FR Us lying below and above, respec- 
tively, the fiducial luminosity PmMHz ~ 10 25 /(#o/70) 2 
W Hz -1 . At higher radio frequencies the division is 
at Pi.4GHz ~ 2 x 10 24 /(iJ /70) 2 W H z" 1 , with some 
depen dency also on optical luminosity (|Owen fc White! 
1991) and therefore a somewhat large overlap. Finally, 
SFG can also be relatively strong radio emitters but are 
hosted by spiral and irregular galaxies at the relative low 
rcdshifts of our sample (see §[272]). These star forming 
radio sources dominate the local (z < 0.4) radio lumi- 
nosity function below Pi 4GHz ~ 10 23 W Hz" 1 but reach 
only Pi.4qHz ~ 10 24 W Hz" 1 (e.g.. ISadler elaTI [200l 
iBest et al.ll2005D . as compared to Pi. 4 gh z « 10 26 W Hz" 1 
or more for the powerful radio-galaxies. In terms of the 
rest-frame ratio of radio to optical flux density, which we 
define as R = log(Si. 4 GHz/Sv) ( wnere Sy is the V-band 
flux density), SFG have relatively small values. For ex- 
ample , the major i ty of the ultra-luminous infrared galax- 
ies m IVega et al.1 (l2008h have Si. 4GHz ~ Si um , which cor - 
responds to R « 0.4, while iMachalski fc Condon! (|1999l ) 
have found that 98% of their normal and starburst galax- 
ies have R < 1.3 (converting from their definition). It 
is important to note that these two studies reach only 
z ~ 0.1 and ~ 0.2 respectively, and therefore one can- 
not exclude the possibility that higher redshift sources 
might have larger R values (see § 12. 2[) . We assume for 
the time being R = 1.4 as the maximum fi ducial value for 
SFG. We note that ISevmour et all ([20081 have used the 
(observed) Si.4ghz/S2.2/j ratio as a discriminant between 
AGN and SFG. Although we have K band photometry 
for less than half of our sample, we find a very strong cor- 
relation between observed Si.4GHz/S2.2/i and rest-frame 
Si.4GHz/Sy, from which we infer that Seymour et al.'s 
dividing line corresponds to R = 1.4 ± 0.1, consistent 
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with our own (but see § [4]) . 

About ~ 10% of AGN have, for the same optical power, 
radio powers 3 — 4 orders of magnitude higher than the 
rest. Two different definitions of ra dio-loudness have 
been used in the literat ure (see, e.g., iKellermann et al.1 
ll989HStocke et al.lll992h . The first is based on the rest- 
frame ratio of the (5 GHz) radio to optical flux density, 
with radio- loud sources having values > 10. Assuming 
a r = 0.7 (see below), and converting to 1.4 GHz, the 
commonly accepted dividing line between radio-loud and 
radio-quiet AGN translates then to R ~ 1.4, which is the 
same as the maximum value for SFG adopted above. The 
other definition is based on gigahertz radio power and 
appears to be redshift-dependent, with a dividing line 
ranging from w 5 x 10 23 W Hz -1 for the Palomar Green 
sample (optically bright and at relatively low-redshift) 
up to 5 x 10 25 W Hz -1 for the Large Bright Quasar Sur- 
vey sample (optically fainter and at higher redshift ) . The 
two definitions overlap somewhat (|Padovanil [l993). The 
radio power limit is, not surprisingly, very close to the 
maximum power which can be reached by SFG. In the 
following we will use a single value, Pi.4ghz ~ 10 24 W 
Hz -1 , appropriate for our relatively low redshift sample 
(see below), and comment also on how our results depend 
on this value. 

2.2. The CDFS-VLA sample 

Our sample includes all CDFS-VLA sources with a 
highly reliable (likelihood ratio > 0.2 as defined in Pa- 
per II) optical counterpart (apart from one source very 
close to a bright star for which we could not get reli- 
able photometry) and eight empty fields, for a total of 
256 objects, 194 of which belong to the complete sample. 
We include all sources when studying the overall sam- 
ple properties but we obviously restrict ourselves to the 
complete sample when dealing with number counts and 
fractions of sources of a given class. 

Our usage of V band data to derive the radio-to-optical 
flux density ratio is dictated by two facts: 1. almost 
complete coverage for our sources (we estimated V mag- 
nitudes from other bands for 20 sources using mean col- 
ors); 2. consistency with previous definitions of radio- 
loud AGN. However, our results are broadly independent 
of this choice. 

One hundred and eight (42%) sources have spectro- 
scopic redshifts, while photometric redshifts are also 
available for 78 (31%) objects (see details in Paper II). 
In total, we have redshift information for 73% of our ob- 
jects, with this fraction reaching 77% for the complete 
sample. The distribution peaks at z ~ 0.7, with a sharp 
cut-off at z > 1 (see Fig. 11 in Paper II). The mean 
redshift is (z) = 0.80 ± 0.05 (median 0.67). 

Fig. Q]shows the V band magnitude versus the 1.4 GHz 
radio flux density for the CDFS-VLA sample. The lines 
represent different values of the logarithm of the radio- 
to-optical flux density ratio, ranging from 5 (top right) 
down to —1 (bottom left). We note that 63% of the 
sources have observed R values beyond the star forming 
regime, with some of the sources reaching extremely high 
values of R, equivalent to radio-to-optical flux density 
ratios w 10 5 . 

This is better seen in Fig. [2] (top panel) , which shows 
the distribution of rest-frame R. For K-correction pur- 
poses we used the observed radio spectral index, when 
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Fig. 1. — The V band magnitude versus the 1.4 GHz radio flux 
density for the CDFS-VLA sample. Arrows represent lower limits 
of V m ag • The lines represent different values of the logarithm of the 
observed radio-to-optical flux density ratio R = log(Si.4GHz/'Sv), 
ranging from 5 (top right) down to — 1 (bottom left). The dashed 
line indicates the maximum local value for SFG and the ap- 
proximate dividing line between radio-loud and radio-quiet AGN 
(Jl = 1.4). 

available, and the mean value, (a r ) ~ 0.7, when that 
was not the case. In the optical band K-correction val- 
ues appropriate for the morphological type of t he objects 
were adopted, based on IColeman et al.1 (|1980f ). Sources 
without morphological classification were assumed to be 
ellipticals 1 (see below), while for compact objects (likely 
AGN), a = 0.5 was assumed. The mean value is 
(R) = 1.04 ± 0.06. The distribution of R for the objects 
without redshift (derived assuming z = (z)) is widely dif- 
ferent (P > 99.99% according to a Kolmogorov-Smirnov 
test) from that of the remaining sources. The mean val- 
ues for the two sub-samples are in fact (R) — 1.62 ± 0.12 
and (R) = 0.83 ± 0.06 respectively, with 60% of the 
sources without redshift being in the star forming regime. 
Their larger R values are due to their fainter V magni- 
tudes, as the radio flux density distributions for the two 
sub-samples are not significantly different. 

Fig. [2] (lower panels) shows the R distribution for vari- 
ous morphological classes (see Paper II for more details), 
namely ellipticals and lenticulars (51), spirals (60), ir- 
regular galaxies (including mergers; 42), compact ob- 
jects (8), and sources which are either not easily clas- 
sifiable, too faint, or complex (95). The various types 
have somewhat different R values, as shown in Table [TJ 
Compact sources have the largest (R) value, followed by 
irregulars, and then galaxies which are not easily clas- 
sifiable. Ellipticals/lenticulars and spirals have the low- 
est radio-to-optical flux density ratios. The R distribu- 
tion for galaxies not classifiable is significantly different 

1 Note that, given the larger K-corrections of ellipticals, any 
other assumption would result in larger R values. 
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Fig. 2. — The distribution of the rest-frame radio-to-optical flux 
density ratio R = log(Si.4GHz/'S'v)- Top panel: Full sample. The 
dashed area represents lower limits. The vertical dashed line indi- 
cates the approximate maximum value for local SFG and dividing 
line between radio-loud and radio-quiet AGN. Other panels: Ob- 
jects sorted by morphological classes, namely, from top to bottom: 
ellipticals and lcnticulars, spirals, irregulars, compact, and sources 
which are either not easily classifiable, faint, or complex. 

(P ~ 98.5%) from that of spiral/irregular galaxies and 
also from that of cllipticals/lenticulars (P > 99.99%). 
Somewhat surprisingly, spirals and irregulars, which we 
can group together as likely hosts of SFG, have R larger 
than ellipticals/lcnticulars, with the two R distributions 
being significantly different (P ~ 98.8%). Only 12% of 
elliptical/lenticular galaxies have R > 1.4, while 21% of 
spirals and irregulars are above this value. We find a 
significant difference also between spirals and irregulars 
(P > 99.99%). It has been suggested that the radio-to- 
optical flux density ratio alone canno t be used to discrim- 
inate between AG N and SFG (e.g.. lAfonso et al.ll2005l 
iBarger et al.l 12007). We confirm that this is indeed the 
case. 

Possible causes for the fact that spirals and irregu- 
lars in our sample reach relatively high radio-to-optical 
flux density ratios, include: a) R evolution. Our sam- 
ple has (z) ~ 0.8, reaching z ~ 2, as compared to the 
local star forming samples for which R has been pre- 
viously derived. We have quantified possible evolution- 
ary effects by assuming a radio evolution cx (1 + z) 2 7 
(Hopkins 2004 ) and passive V-ban d evolution for Sa and 
Sc galaxies (jColeman et al.lll980t ). In the Sa case R is 
predicted to decrease at higher redshifts, because the op- 
tical evolution is stronger than the radio one. For Sc 
galaxies, on the other hand, there is a weak evolution 
with R{z) oc log(l + z). This means that at z = (z) 
-ffmax < 1-7 and at z ~ 2 i? max should be < 1.9. This ef- 
fect could explain up to half of the R > 1.4 sources with 
redshift; b) uncertainties in morphological classification. 
While Paper II has shown that overall our morphological 
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Fig. 3. — The rest-frame radio-to-optical flux density ratio R = 
log(S'i AGHz/ Sv) versus the 1.4 GHz radio power for the sources 
with redshift information. The horizontal (long-dashed) line at 
R = 1.4 and the vertical (short-dashed) line at log P T = 24 indicate 
the approximate maximum values for local SFG and dividing lines 
between radio-loud and radio-quiet AGN respectively. Sources are 
classified morphologically, namely: ellipticals and lenticulars, red 
filled circles; compact, blue triangles; spirals, green filled squares; 
irregulars, cyan open squares; sources which are either galaxies not 
easily classifiable, faint, or complex, black open circles. The black 
open squares with error brackets represent the locus of objects 
without redshift assuming (from left to right) z = (z), z = 1.0, 
z = 1.5, and z = 2 and an optical K-correction appropriate for 
ellipticals. See text for more details. 

classification is robust, there could still be subtle effects, 
which might introduce some bias. For example, while 
the distributions of magnitudes of ellipticals/lenticulars 
and spiral galaxies in our sample are very similar, irreg- 
ulars are on average two magnitudes fainter. If some 
galaxies are preferentially classified as such when they 
are faint, this will result in artificially large R values for 
irregulars, with the effect being amplified by the fact that 
the K-correction for irregulars is the smallest one; c) K- 
correction. We assumed standard K-corrections, which 
imply that all spectral energy distributions (SEDs) for 
sources of a given class are the same. This is obviously 
unrealistic and will certainly introduce a scatter in the 
R distributio n; d) the pre s ence of sources like the one 
discovered by iNorris et alj (|2008t ): a classic triple radio 
source at z = 0.93 and P r ~ 4 x 10 25 W Hz" 1 with the 
SED of a spiral galaxy. Finally, it is important to notice 
that not all spirals and irregulars are SFG, as potentially 
many of these sources could be AGN. It is then clear that 
we need to add luminosity information. 

Fig. El shows R versus the 1.4 GHz radio power for 
the sources with redshift information. Apart from the 
obvious correlation between R and P r (the former be- 
ing equal to P r /P opt ), Fig. O shows that all but three 
sources with R > 1.4 have P T > 10 23 W Hz" 1 . Turn- 
ing to morphological source classification, the following 
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TABLE 1 

Mean values. 



Galaxy class 




(R) 


N 


{logPr-} 


N 


(log L x ) 


Ellipticals and lenticulars 


0.52 ±0.11 


51 


23.47 ±0.16 


47 


41.55 ± 0.23 


Spirals 




0.58 ± 0.09 


60 


22.75 ±0.12 


59 


40.65 ±0.16 


Irregulars 




1.42 ±0.09 


42 


23.38 ±0.11 


36 


41.49 ± 0.27 


Spirals and Irref 


;ulars 


0.92 ± 0.07 


102 


22.99 ±0.09 


95 


40.87 ±0.15 


Compact 




2.40 ±0.27 


8 


24.62 ±0.66 


6 


43.54 ±0.55 


Other 




1.34 ±0.10 


95 


24.05 ±0.13 


38 


42.55 ± 0.16 


No redshift, z = 


(*) 


1.62 ±0.12 


70 


23.92 ±0.09 


70 


41.85 ± 0.13 


No redshift, z = 


1.0 


1.31 ±0.12 


70 


24.15 ±0.09 


70 


42.06 ± 0.13 


No redshift, z = 


1.5 


0.83 ±0.13 


70 


24.56 ± 0.09 


70 


42.41 ± 0.13 


No redshift, z = 


2.0 


0.68 ±0.14 


70 


24.85 ±0.09 


70 


42.68 ±0.13 



points can be made: a) the large majority (85%) of spi- 
rals and irregulars have P r < 10 24 W Hz -1 , with this 
fraction reaching 97% and 100% for P r < 3 x 10 24 W 
Hz -1 and P r < 10 25 W Hz -1 respectively; b) spirals and 
irregulars with R > 1.4 have relatively high radio pow- 
ers ((P r ) - 6 x 10 23 W Hz" 1 ) and basically all of the 
sources with R > 2 are above (or very close to) the 10 24 
W Hz -1 limit; c) 5/6 of compact sources are consistent 
with being radio-loud AGN; d) ellipticals and lenticu- 
lars span a very large range of radio powers, consistent 
with that of radio-galaxies; there is substantial overlap 
at low powers between early- and late-type galaxies, with 
75% of the former having P r < 10 24 W Hz" 1 ; e) all 
"other" sources have powers in the range of radio galax- 
ies and half of them actually reach into the radio-loud 
AGN regime. Most of the various types have quite dis- 
tinct radio powers, as shown in Table[TJ Compact sources 
have the largest radio powers, followed by galaxies which 
are not easily classifiable, and ellipticals/lenticulars. Spi- 
rals/irregulars have the lowest values. The radio power 
distribution for the sources without morphological classi- 
fication is significantly different (P > 99.99%) from that 
of spirals and irregulars and also (P ~ 98.8%) from that 
of early- type galaxies. 

We have also o ptical spectral classification, from 
ISzokolv et all (|2004f ). for a relatively small fraction (14%) 
of the sample, which we use to gather further informa- 
tion. This includes Broad Line AGN (BLAGN), high ex- 
citation line objects (HEX), low excitation line objects 
(LEX), Narrow Emission Line Galaxies (NELG), and ob- 
jects with normal galaxy spectra (ABS). Broadly speak- 
ing, BLAGN correspond to type-1 AGN, HEX mostly to 
type-2 AGN, while the properties of NELG, LEX, and 
ABS are consistent with those of normal galaxies but in 
many cases this is simply a limitation of the optical data. 

The spectral classification adds the following informa- 
tion: a) AGN have, on average, radio powers larger than 
those of LEX, ABS, and NELG (~ 10 24 W Hz" 1 as com- 
pared to - 8 x 10 22 W Hz" 1 ); b) all but one LEX, ABS, 
and NELG have P r < 10 24 W Hz" 1 and - 78% of them 
concentrate in the bottom-left quadrant of Fig. [3l con- 
sistent with the fact that this is where SFG are located. 

As for the radio-to-optical flux density ratio, evolution 
will also affect t he maximum r adio power of SFG. For 
P T oc (1 + z) 2 - 7 (|Hopkindl200l . this could reach P r ~ 
5 x 10 24 W Hz" 1 and P r ~ 2 x 10 25 W Hz" 1 at z = (z) 
and z ~ 2 respectively. 

Fig. [3] does not include sources without redshift. If we 



assign a redshift equal to the mean value of the sam- 
ple to such sources, they would have (R) ~ 1.6 and 
logP r ~ 23.9 (Table [TJ, that is most of them would 
fall into a region outside that of SFG (see Fig. [3|). For 
z = 1 — 2, which are more likely values given the de- 
pendence we find between redshift and magnitude, radio 
powers would obviously increase, with R values decreas- 
ing due to the larger K-correction (as ~ 80% of such 
sources have no morphological information and we as- 
sume that they are ellipticals), as shown in Fig. [31 How- 
ever, the point remains that the majority of these sources 
would not have values very typical of SFG. 

In summary, R is not a very good discriminant be- 
tween SFG and AGN. Radio power fares somewhat bet- 
ter but the overlap between the two classes is still quite 
large. This is obviously particularly bad for the sources 
for which no optical morphology is available. 

3. ADDING X-RAY INFORMATION 

We now consider the X-ray information to learn more 
about the nature of the radio sources. Fig. [4] plots the 
rest-frame radio-to-optical flux density ratio R vs. the 
2 — 10 kcV power for the sources with redshift informa- 
tion. X-ray powers and upper limits are available for 
all but three sources and for all of those in the com- 
plete sample (see Paper III). The vertical short-dashed 
line at L x = 10 42 ergs s _1 indicates the approximate 
dividing line betwe en AGN and normal galaxies (e.g., 
ISzokolv et aT][200l . 

To properly take into acco unt the upper limits on X- 
ray power we used ASURV (|Lavallev et al.| [l992). the 
Survival An alysis package which emp loys t he routines 
descri bed in iFeigelson fc Nelson! (|1985h and llsobe et al.l 
(1986), which evaluate mean values by dealing properly 
with non-detections. Survival analysis methods also eval- 
uate probability of correlation and linear regression fits 
by dealing properly with non-detections. A correlation 
found between luminosities, when such non-detections 
are included, is a true correlation, while flux density - 
flux density diagrams, which are sometimes used, will re- 
veal the intrinsic relationship between luminosities only 
under very limited circumstances (no censored points, 
linear relationship). Since upper limits make up the ma- 
jority (~ 67%) of the sample, the uncertainties on our 
results are likely to be larger than the ones given by the 
software package. 

We find a very strong (P > 99.0%) correlation between 
I/ x and R for the whole sample, which is even stronger 
(P > 99.9%) for spirals and irregulars. Two striking 
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Fig. 4. — The rest-frame radio-to-optical flux density ratio R = 
log(S 1 4QH z /5y) versus the 2 — 10 keV X-ray power for the sources 
with redshift information. The horizontal (long-dashed) line at 
R = 1.4 indicates the approximate maximum value for local SFG 
and dividing line between radio-loud and radio-quiet AGN, while 
the vertical (short-dashed) line at L x = 10 42 ergs s — 1 indicates 
the dividing line between AGN and normal galaxies. Sources are 
classified morphologically, namely: ellipticals and lenticulars, red 
filled circles; compact, blue triangles; spirals, green filled squares; 
irregulars, cyan open squares; sources which are either galaxies not 
easily classifiable, faint, or complex, black open circles. The black 
open squares with error brackets represent the locus of objects 
without redshift assuming (from left to right) z = (z), z = 1.0, z = 
1.5, and 2 = 2 and 2 = 2 and an optical K-correction appropriate 
for ellipticals. See text for more details. 



features of Fig. H] are: a) the fact that 58% of all X- 
ray detections have L x > 10 42 ergs s _1 , with this frac- 
tion reaching 62% if one includes upper limits as well; 
b) the almost complete lack of objects in the top left 
quadrant, as only 3% of the sources have R > 1.4 and 
L x < 10 42 ergs s _1 (only half of those are irregulars). 
Indeed, R > 1.4 sources have (logi x ) = 42.14 ±0.28, al- 
most an order of magnitude larger than R < 1.4 sources, 
which have (logL x ) =41.21 ±0.13. This suggests that 
a large fraction of our sources are AGN and that many 
high R sources with spiral/irregular morphology cannot 
be SFG as their X-ray powers are too high. 

As was the case for the radio power, most of the var- 
ious types have quite distinct X-ray powers, as shown 
in Table [T] Compact sources have the largest X-ray 
powers, followed closely by galaxies which are not easily 
classifiable, and ellipticals/lenticulars. Irregulars/spirals 
have the lowest values. The X-ray power distribution 
for spirals and irregulars is significantly different 2 from 
that of the sources without morphological classification 
(P > 99.98%) and also from that of early-type galaxies 
(P > 98.6%). The spectral classification adds the fact 

2 The ASURV package evaluates the difference between two dis- 
tribution using five different tests. We quote here the lowest prob- 
ability. 



that AGN, on average, have X-ray powers larger than 
those of LEX, ABS, and NELG (~ 2 x 10 43 ergs s" 1 as 
compared to 6 x 10 41 ergs s -1 ), although ~ 30% of LEX 
and ABS sources are in the AGN region (all BLAGN 
and all but but two HEX have AGN-like X-ray powers); 
in particular, all NELG have only upper limits on their 
X-ray fluxes. 

As before, Fig. [4] does not include sources without 
redshift. If we assign a redshift equal to the mean value 
of the sample to such sources, they would have (R) ~ 1.6 
and logi x ~ 41.8 (Table[T]), that is many of them would 
fall into a region outside that of SFG. For z = 1 — 2, 
which are more likely values given the dependence we 
find between redshift and source magnitude, R would 
decrease but X-ray powers would obviously increase, as 
shown in Fig. 01 with the result that the majority of 
these sources would not have values typical of SFG. 

Fig. M plots 2 — 10 KeV X-ray power versus the 
1.4 GHz radio power for the sources with redshift in- 
formation. The horizontal, long-dashed line at L x — 
10 42 ergs s _1 indicates the approximate dividing line 
betw een AGN and normal galaxies (e.g., ISzokolv et al.l 
2004), while the vertical, short-dashed line at logP r = 
24 indicates the approximate maximum radio power 
(at 2^0) due to star-formation, which roughly co- 
incides with the dividing lines between low-redshift 
radio-loud and radio-quiet AGN. The dot-long-dashed 
line is the X-ray-radio power re lation for nearby SFG 
(|Ranalli. Comastri fe Settil [2003). the dotted lines rep- 
resent the approxi mate locus of radio galaxie s in th e 
AGN catalogue of iPadovani, Giommi fe Fiora (|1997h . 
while the dot-short-dashed and short-long-dashed lines 
are t he best fits for radio-qu iet and radio-loud quasars 
from iBrinkmann et al.l (|2000T ). where their 2 keV pow- 
ers were converted to the 2 — 10 keV band assuming 
a x = 0.8. For K-correction purposes we used the ob- 
served radio and X-ray spectral indices, when available, 
and the mean values, (a r ) ~ 0.7 and (a x ) ~ 0.8, when 
that was not the case. 

Fig. [5] shows that there is a correlation between ra- 
dio and X-ray luminosity for the whole sample. ASURV 
confirms this with a very high significance (P > 99.99%) 
giving L x oc p0- 93±0 - 08 . Basically all sources in the radio- 
loud AGN region have AGN-like X-ray powers. The 
opposite, however, is not true, that is not all sources 
with L K > 10 42 ergs s" 1 have P r > 10 24 W Hz" 1 . This 
makes sense, as not all AGN are radio-loud. Turning to 
morphological source classification, the following points 
can be made: a) 38% of spiral and irregular galaxies 
have X-ray detections inconsistent with the powers typ- 
ical of galaxies (> 10 42 ergs s _1 ); that is, a sizable frac- 
tion of these sources are most likely powered by AGN. 
The best fit correlation using ASURV (green line) ex- 
cludin g sources with P r > 10 23 5 WHz -1 for consistency 
with iRanalli. Comastri fc Settil (|2003[ ) and X-ray detec- 
tions with L x > 10 42 ergs s _1 to select against AGN, is 
L x oc p0- 8±0 1 (P > 99.6%). Given the relatively large 
fraction of X-ray upper limits (~ 70%), we regard this 
as not inconsistent with the X-ray-radio power relation- 
ship for nearby SFG (dot-dashed line); b) as was the 
case for the radio power, ellipticals and lenticulars span a 
very large range of X-ray powers, consistent with that of 
radio-galaxies (red line) ; again, there is substantial over- 
lap at low X-ray powers between ellipticals and spirals; 
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Fig. 5.— The 2-10 keV X-ray power vs. the 1.4 GHz ra- 
dio power for sources with redshift and X-ray information. The 
horizontal (long-dashed) line at L x = 10 42 ergs s — 1 indicates 
the dividing line between AGN and normal galaxies, while the 
vertical (short-dashed) line at P r = 10 24 W Hz -1 shows the 
maximum P T of SFG at z ~ 0, which coincides with the di- 
viding lines between low-2 radio-loud and radio-quiet AGN. The 
dot-long-dashed line is the L x — P T relationship for nearby SFG 
(Ranalli, Comastri & Sctti 2003), converted to our cosmology, the 
dotted li nes denote the locus of rad i o gala xies in the AGN cat- 
alogue of Padovani. Giommi & Fiorc (1997), while the dot-short- 
dashed and short-long-dashed lines are the best fi ts for radio-quiet 
and radio-loud quasars from Brinkmann ct al. (2000). Arrows 
represent X-ray upper limits. Sources are classified morphologi- 
cally, namely: ellipticals and lenticulars, red filled circles; com- 
pact, blue triangles; spirals, green filled squares; irregulars, cyan 
open squares; sources which are either galaxies not easily classi- 
fiable, faint, or complex, black open circles. Solid lines show the 
best fits in each morphological category taking the upper limits 
into account (up to logP r < 23.5 and excluding X-ray detections 
with logL x > 42 for spirals/irregulars: green line). 

c) "other" sources have mostly L x > 10 42 ergs s _1 , with 
the best fit correlation (black line) into the radio-galaxies 
area; d) compact sources have X-ray powers > 5 x 10 41 
ergs s _1 and » 10 times larger than the other classes at 
a given radio power. 

Note that th e correlation we find for spirals/irregulars 
disagrees with iBarger et al.l (|2007l) , who suggested that 
the Ranalli et al. correlation was due to selection effects, 
but agrees with other results (see discussion in Paper 
III). 

We also notice that only a relatively small number of 
sources cluster around the relationships previously de- 
rived for radio-quiet and radio-loud quasars. We will 
return to this point later in § [4] Here we point out that 
the fraction of QSOs, defined as sources with L x > 10 44 
ergs s _1 , is only 7% of the objects with redshift and 
X-ray information. Note that this has to be consid- 
ered a robust upper limit, as radio galaxies are also 
strong X-ray emitters. For example, the relatively low- 
redshift ( (z) = 0.1) radio galaxies i n the AGN cat- 
alogue of iPadovani, Giommi fc Fiorel (|1997[ ) can reach 



L x ~ 3 x 10 44 ergs s 1 , and all fiv e 3CR radio galax- 
ies studied by iSalvati et all (|2008t ) ((z) = 1.3) have 
L x > 10 44 ergs s _1 . BL Lacs are even stronger X-ray 
e mitters, with ~ 70% of tho s e in t he AGN catalogue 
of IPadovani. Giommi fc Fiorel (|1997h having L x > 10 44 
ergs s _1 , but are likely to be a minority population in 
our sample (see § I4.3j) . Indeed, the average radio spec- 
tral index for the L x > 10 44 ergs s _1 sources is ~ 0.6, 
with only ~ 20% of them having a T < 0.5 and therefore 
qualifying as flat-spectrum, which is one of the defining 
properties of BL Lacs. 

In summary, the X-ray data appear to be quite power- 
ful in recognizing AGN among the sub-mJy population. 

4. THE SUB-MJY SOURCE POPULATION 

4.1. The redshift sub-sample 

We now examine collectively the radio, optical, and 
X-ray data. We start by considering only sources with 
redshift information, which include 77% of the complete 
sample. Based on the results presented in the previous 
sections, our candidate star-forming galaxies are defined 
as fulfilling all the following requirements: 

1. R< 1.7 

2. P T < 10 24 - 5 W Hz" 1 

3. optical morphology different from elliptical or 
lenticular 

4. L x < 10 42 ergs s _1 for X-ray detections, no limit 
otherwise. 

The first two criteria include 87% of spirals and irregu- 
lars, with the remaining fraction being most likely AGN, 
the third one excludes sources not associated with star 
formation at the redshifts under consideration, while the 
fourth one excludes AGN. 

Objects not fulfilling these criteria were considered to 
be AGN. Note that the last requirement is conservative 
in the sense that it maximizes the number of SFG, as 
some of the sources with X-ray upper limits above the 
X-ray threshold could still have L x > 10 42 ergs s _1 (see 
below). 

Radio-quiet AGN were also selected by considering 
sources with R < 1.4 and X-ray detections with L x > 
10 ergs s _1 . We did not apply any cut in radio power 
since the objects with P r > 10 24 W Hz -1 are all at rela- 
tively high redshifts ((z) ~ 2) and are therefore expected 
to have relatively large luminosities even if radio-quiet 
(see §[23]). 

Table [2] and Fig. [6] present the Euclidean normalized 
number counts for the total redshift sub-sample and also 
for SFG and AGN, derived as described in Paper I. AGN 
make up 4ltg% of sub-mJy sources with redshift infor- 
mation but their number drops towards lower flux den- 
sities, going from 100% at ~ 10 mJy down to 34% at 
the survey limit. SFG, on the other hand, constitute 
59lg°% of the sub-mJy sample having measured red- 
shifts. SFG are missing at high flux densities but be- 
come the dominant population below w 0.1 mJy, reach- 
ing 66% at the survey limit. Radio-quiet AGN represent 
23± 6 5 % (or - 56% of all AGN) of sub-mJy sources but 
their fraction appears to increase at lower flux densities, 
where they make up 81% of all AGN at the survey limit. 
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TABLE 2 



Euclidean normalized 1.4 GHz counts for sources with redshift. 
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Fig. 6. — The Euclidean normalized 1.4 GHz CDFS source counts 
for sources with redshift information: total counts (black triangles), 
star forming galaxies (filled green circles), all AGN (red squares), 
and radio- quiet AGN (op en blue circles). Error bars correspond to 
la errors ( Gchrcls 1982). 



TABLE 3 

Fraction of sub-mJy star forming galaxies for various 
parameter limits. 
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log L x 
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We caution that the real uncertainties on the radio- 
quiet AGN counts are larger than indicated owing to 
two effects, which go in opposite directions: 1) ~ 60% 
of our putative radio-quiet AGN have elliptical or unde- 
fined morphology, which means their optical fluxes have 
been K-corrected assuming an elliptical spectrum. This 
could result in an artificial lowering of R if the AGN 
component is non negligible in the optical band, trans- 
lating in turn in an overestimation of their number (as 
we selected sources with R < 1.4). Our selection crite- 
ria could also include some radio galaxies, although their 
number is expected to be small; 2) the inclusion of X- 
ray detections only could underestimate the number of 
radio-quiet AGN, as some sources with upper limits on 
their X-ray power above 10 42 ergs s -1 could still have 
intrinsic powers above this limit. 

To assess the dependence of our results on our selec- 
tion criteria we show in Table [3] how the fraction of sub- 
mJy SFG changes with different parameters. Namely, 
we changed the values of R (by +0.6, —0.3), P r (by ±0.5 
dex), and L x (by ±0.5 dex). The fraction of sub-mJy 
SFG changes remarkably little, with a maximum increase 
and decrease of +3% {R = 2.3) and -6% (i? = 1.4) re- 
spectively, well within the la error for our standard set 
of values. The biggest change happens when we include 
in our selection only sources with upper limits on their 
X-ray powers < 10 42 5 ergs s _1 , instead of all the X-ray 
non detections, to check for the effect of our assumption 
that all sources with X-ray upper limits fall below the 
AGN dividing line. The fraction of SFG decreases by 
almost one third to 41 ± 7% (a 1.6a difference), which 
shows that indeed our assumption is a conservative one. 

4.2. The whole sample 

We now extend our analysis to the whole CDFS-VLA 
sample. As discussed in § 12.21 and § [3] (see also Figs. [3] 
and[4|), in fact, the bulk of the sources without redshift 
are unlikely to be SFG. Therefore, it would be unwise to 
simply extrapolate the results of the redshift sub-sample 
to the whole sample. For the same reason, any compari- 
son with predicted number counts needs to be done with 
the full sample. 

We then apply our selection criteria described in § 14.11 
to all our sources, making various assumptions on the 
missing redshifts. Table 0] shows that the fraction of 
sub-mJy SFG is only weakly dependent on the assumed 
redshift for the sources without spectroscopic and pho- 
tometric redshifts. This is due to the fact that the in- 
crease in power for larger redshifts is compensated by 
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TABLE 4 
Fraction of sub-mJy star 
forming galaxies for 
different assumed redshifts. 



rcdshift 


N 


sub-mJy fraction 
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88 
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58 ±8 


2.0 
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the decrease in R due to the larger K-correction, since 
we assumed elliptical morphology for the sources without 
that information. For z = (z) the fraction of sub-mJy 
SFG equals 57 ± 8% and in the following we make this 
assumption. 

Table [5] and Fig. [7] present the Euclidean normalized 
number counts for the whole sample and also for SFG 
and AGN, derived as described in Paper I. The trends 
are very similar to those of the redshift sub-sample. AGN 
make up 41 ± 6% of sub-mJy sources and their counts are 
seen to drop at lower flux densities, going from 100% of 
the total at ~ 10 mJy down to 36% at the survey limit. 
SFG, on the other hand, which represent 57 ± 8% of the 
sample, are missing at high flux densities but become the 
dominant population below « 0.1 mJy, reaching 64% at 
the survey limit. Radio-quiet AGN represent 19l 4 % (or 
46% of all AGN) of sub-mJy sources but their fraction 
appears to increase at lower flux densities, where they 
make up 66% of all AGN and w 24% of all sources at the 
survey limit, up from w 10% at ~ 1 mJy. 

4.3. Expected number of radio sources 

In Fig. [7] we overlay some model calculations for the 
AGN and SFG contribution to the radio source counts. 
The AGN model includes FR I and FR II radio galax- 
ies and BL Lacertae objects, the only populations sig- 
nificantly present in the counts below a few mJy at 1.4 
GHz (Padovani, in preparation). In fact, according to 
unified schemes, radio quasars, being FR lis with their 
powers boosted by relativistic beaming, are expected to 
disappear from the radi o population below S T < 1 mJy 
(jPadovani fc UrrvlflQOl . 

Radio galaxy counts were calculated starting from the 
2.7 GHz luminosity functions of 2 Jy FR I a nd II radio 
galaxies derived by lUrrv fc Padovanil (|1995l ) , converted 
to 1.4 GHz assuming a r — 0.7. Based on their results, 
no evolution was assumed for FR Is, while their best fit 
luminosity evolution (t = 0.26 for an exponential evo- 
lution P(z) = P(0)exp[T(z)/r], where T(z) is the look- 
back time and r is the evolutionary time scale in units 
of the Hubble time) was used for FR lis. BL Lac num- 
ber counts were estima ted from the beaming model of 
lUrrv fc Padova ni ( 1995) (converted from 5 GHz assum- 
ing a r = —0.1), which fits very well not only the original 
1 Jy data on which it w as based but also the de eper (~ 50 
mJy) DXRBS sample (|Padovani et alJl2007af l. 

The radio AGN counts are dominated by FR I radio 
galaxies at sub-mJy levels. In fact, given their higher 
radio power and stronger evolution, FR lis reach only 
S r rs 0.5 mJy. FR Is appear to make up ~ 40% of 
the radio galaxies with S T >1 mJy but ~ 90% of the 
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Fig. 7. — The Euclidean normalized 1.4 GHz CDFS source 
counts: total counts (black triangles), SFG (filled green circles), 
all AGN (red squares), and radio-q uiet AGN (ope n blue circles). 
Error bars correspond to la errors {Gchrcls 1986). Model calcu- 
lations refer to SFG (green dot-dashed lines), displayed with a lo" 
range on the evolutionary parameters, all AGN (red dashed line), 
radio-quiet AGN (blue dotted line), and the sum of the first two 
(black solid line). See text for more details. 



whole sample, and to reach flux densities as low as ~ 
0.1 /xJy. BL Lacs can also reach quite faint flux densities 
(i~ 10 fiJy) but make up only a minority of the AGN, 
« 10 — 15% in the range 0.1 — 1 mJy. 

The su r face d ensity for SFG has been derived following 
iHopkinsI (l2004h. Namely , the local luminosity function 
from iSadler et all lj2002T) was adopted, with luminosity 
evolution oc (1 + z) 2 ' 7 to z — 2 (and constant thereafter) 
and density evolution oc (1 + z) 015 . To have a feeling 
for the likely uncertainties, we have also evaluated the 
number counts by varying the evolutionary parameters 
by la. 

Our AGN counts include a radio-quiet component, 
as radio-quiet AGN are not radio-silent and are rou- 
tinely detected at faint flux densities in the radio band. 
To estimate their expected number counts one would 
need a sizable sample of such radio-selected sources from 
which to derive the radio luminosity function and evo- 
lution, to then predict their number densities at lower 
flux densities, as we have done for radio- and star form- 
ing galaxies. Unfortunately, such a sample is not avail- 
able. Radio-quiet AGN, however, are very strong X- 
ray emitters and indeed dominate the X-ray sky. We 
th erefore took the 5 — 10 keV AGN nu mber counts give n 
bv lComastri fc Brusal (|2008h (based on lGilli et aill2007h . 
which fit the observed X-ray counts down to ~ 10 -15 ergs 
cm~ 2 s _1 and make predictions down to ~ 10 -17 ergs 
cm~ 2 s _1 , to obtain an estimate of the radio-quiet AGN 
counts in the radio band. We first multiplied the X-ray 
counts by 0.9, to take into account the fact that radio- 
loud AGN make up « 10% of all AGN, and then con- 
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TABLE 5 

Euclidean normalized 1.4 GHz counts for the whole sample. 
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verted them to Euclidean normalized counts. To derive 
radio number counts one would then need to know the 
intrinsic radio-to-Xray flux ratio distribution, or at least 
the mean ratio for radio-quiet AGN. This is at present 
unknown. However, one can obtain robust bounds on 
this ratio by looking at X-ray selected samples and our 
o wn sources. 

iPolletta et all (|2007l ) present the SEDs of a hard X- 
ray selected sample of AGN with /2-iokcV > 10~ 14 ergs 
cm~ 2 s _1 . These include radio flux densities, mostly 
upper limits. We used survival analysis to derive a k- 
corrected mean value {Si.^GHz / h-iokev) = 14.57 ±0.20 
(where the radio flux density is in fiJy and the X-ray 
flux is in c.g.s. units). Radio- loud AGN, which we de- 
fine as sources with L r > 10 24 W Hz -1 and make up 
~ 13% of the sample, were excluded. The radio-quiet 
AGN in our complete sample, defined as done in § 14.11 
have (<5i.4GHz//2-iokev) = 16.72 ± 0.16. Both of these 
values are heavily influenced by selection effects. X- 
ray selection will tend to favor sources with relatively 
small radio-to-X-ray flux ratios, while the opposite will 
be true for radio selection. The truth, as always, will be 
somewhere in the middle and actually quite far from the 
two extremes. We then take as our best guess for the 
radio-to-X-ray flux ratio for radio-quiet AGN the mean 
of the two values, that is (Siaghz I /2-iokev) = 15.6, be- 
ing fully aware that this will only give an order of mag- 
nitude estimate. We note that this value is the same 
as tha t obtained from the sample of Brinkmann et alJ 
(2000) by converting their 2 keV powers to the 2—10 
keV band assuming a x = 0.8. Establishing why this is 
so, given the fact that the iBrinkmann etT al. (2000) sam- 
ple was derived from a cross-correlation between ROSAT 
All-Sky Survey and 1.4 GHz FIRST data, is beyond the 
scope of this paper. We just remark that this validates 
the usage of this value as typical in the literature (e.g., 
ISimpson et "all 120061 : IPolletta et all l2007h . Radio num- 
ber counts for radio-quiet AGN were then obtained by 
using the hard X-ray counts and the value derived above 
converted to 5 — 10 keV by assuming a x = 0.8 (i.e., 

{SlAGHz/ j5-10kev) = 16.0). 

The agreement between observed and expected num- 
ber counts for the various sub-classes and for the whole 
sample in Fig. [7] is very good, especially given the many 
orders of magnitude difference between our flux density 
limit and those of the samples from which the luminos- 
ity functions were derived. Even the radio-quiet AGN 
counts show a very good agreement with our observa- 
tional estimates. We point out that the inclusion of a 



distribution in radio-to-X-ray flux ratio would "stretch" 
those counts conserving the total number, with the result 
that the drop at lower flux densities would be somewhat 
mitigated. 

4.4. Heavily absorbed AGN: the meaning of R 

Radio selection gives us the chance of identifing ob- 
scur ed AGN which mig ht escape X-ray detection (see, 
e.g.. iDonlev et aill2005l ). One might think that, since 
their optical flux is heavily absorbed, such sources will 
have quite large radio-to-optical flux density ratios. In- 
deed, it has been suggested that the occurrence of such 
large ratios in micro- Jy sources is indicativ e of a high de- 
gree of extinction (e.g.. iNorris et al.ll2005l ). Radio-quiet, 
obscured AGN (the so-called "type 2" AGN), would 
then appear to have R values very large for their ra- 
dio power and should stick out in a R - radio power 
plot. To test this hypothesis we plot in Fig. [8] mean 
R and logP r values for broad-lined AGN, both radio- 
quiet and radio-loud, Seyfert 2s, and star forming and 
ra dio galaxies. The data come from the AGN catalogue 
of lPadovani. Giommi fc Fiord ( 1997). Given the redshift 
distribution of our sample (§ 12. 2p and to compare sources 
within a similar redshift range, we have restricted our 
selection to z < 1. Individual data from the same cat- 
alogue 3 for ~ 55% of the C o mpto n-thick sources from 
the compilation of IComastril (j2004f ) are also shown. In 
such sources, defined as having hydrogen column val- 
ues Nn > 1.5 
visible above ~- 
called "mildly Compton thick" sources, while for larger 
column densities the entire high-energy X-ray spectrum 
is absorbed. Compton-thick sources are expected to be 
heavily absorbed in the optical band as well, and indeed 
~ 80% of them are classified as Seyfert 2s/QSO 2s or 
radio galaxies. 

An important result is apparent from Fig. [8] X-ray 
and optically absorbed sources fall on the same relation 
as optically unabsorbed ones. In other words, a large 
value of the radio-to-optical flux density ratio does not 
indicate a high degree of extinction but, rather, a large 
radio power. A similar story is told by the Compton- 
thick sources, which fall on the same relation as the 
other sources (see also Fig. |3]). This can be explained 
by the fact that there is a limit to how faint the op- 

3 For consistency with our sources the optical magnitudes used 
in the evaluation of the radio-to-optical flux density ratio refer, as 
much as possible, to the total, and not the nuclear, magnitude. 



10 cm , nuclear radiation is only 
10 keV for N H < 10 25 cm" 2 , the so- 
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Fig. 8. — The radio-to-optical flux density ratio R = 
log(Si .4GHz/^v) versus the 1.4 GHz ra dio power for various classes 
of sou rces, from the AGN catalogue of Padovani, Giommi &: Fiord 
lT997ft : Broad Line AGN (BLAGN), radio-quiet (left-most point) 
and radio-loud (right-most point), filled red circles; Seyfert 2s, filled 
magenta square; SFG, open green circle; and radio galaxies with 
logP r < 25 (left-most point) and logP r > 25 (right-most point), 
open blue squares. Individu al data for ~ 55% of the Compton-thick 
sources in the compilation of Comastri ( 2004) are also shown (black 
triangles). Error bars represent the la dispersion. The horizontal 
(long-dashed) line at R = 1.4 and the vertical (short-dashed) line 
at logP r = 24 indicate the approximate maximum values for local 
SFG and dividing lines between radio-loud and radio-quiet AGN. 
See text for more details. 

tical flux will be, since once the central AGN is com- 
pletely obscured the host galaxy will contribute all of 
the flux. Radio-quiet, broad- lined AGN have R values 
only slightly smaller than Seyfert 2s because, given also 
the redshift cut-off, they are mostly nearby and therefore 
with a strong galaxy component. Radio-quiet quasars at 
all redshifts, which have a stronger AGN-related optical 
component, have (R) ~ 0.54 and logP r ~ 24.1, that is 
below the locus of the other sources. 

We notice that two sources in Fig. [3] (top-left quad- 
rant) appear to be off the main relation by having R val- 
ues larger than expected for their radio powers. These 
are at low redshift ( z ~ 0.06 and 0.18 ). One of them 
has been studied by iTozzi et all (|2006T ) (ID 97 in their 
Table 1), who derived a x = 0.30^;^ and N H ~ 0. The 
other one is a much weaker X-ray source with an X-ray 
hardness ratio, defined as HR = (H — S)/(H + S), where 
H and S are the net counts in the hard and soft band 
respectively, ~ 0.41. This implies Ah > 10 22 cm -2 . In 
both cases, then, there is no strong indications of absorp- 
tion. Both sources have very faint optical magnitudes for 
their redshifts, clearly off the locus of the majority of our 
sources, which explains their high R. This suggests that 
they might be somewhat peculiar. 

Based on all of the above, it appears that radio and 
optical data cannot identify heavily absorbed sources, 
which might have escaped X-ray detection. This can be 



investigated further using Spitzer data. 

5. DISCUSSION 

By using a variety of multi- wavelength information for 
256 radio sources in the CDFS area, we have shed light 
on the long standing problem of the nature of sub-mJy 
radio sources. We find that SFG make up < 60% of such 
sources. We believe our results are quite robust, since 
whenever we made some assumptions, we chose those 
which would maximize the number of SFG. AGN make 
up the remaining fraction (> 40%) of radio sources with 
1.4 GHz flux densities < 1 mJy. The observational evi- 
dence, backed up by model calculations and an extrapo- 
lation from the hard X-ray band, strongly supports the 
suggestion that AGN are roughly equally split into radio- 
loud, mostly low-power radio galaxies, and radio-quiet 
sources. 

The population mix appears to be flux density depen- 
dent, with SFG becoming increasingly relevant at lower 
flux densities until they become the dominant compo- 
nent below as 0.1 mJy, and AGN exhibiting the opposite 
trend. This seems to be driven by the fall of radio-loud 
sources, as radio-quiet ones display relatively flat counts 
(within the rather large errors). SFG do appear to be 
responsible for most of the well-known up-turn in the 
counts, with radio-quiet AGN providing only a relatively 
small (« 25%) contribution. 

The availability of optical morphology and X-ray infor- 
mation (including upper limits) for ~ 70% and all of the 
sources in the complete sample respectively have been 
fundamental to classify our sources. X-ray data, in par- 
ticular, have been vital to eliminate AGN contaminants 
from spirals and irregulars and to classify sources with- 
out morphological information. We have obviated the 
presence of many X-ray upper limits by using them con- 
servatively and thereby providing robust upper limits to 
the fraction of SFG. Deeper X-ray data should allow us 
to better define the fraction of SFG. 

As mentioned in the Introduction, our results differ 
somewhat with the many papers, which over the years 
have suggested a large dominance of SFG at sub-mJy 
levels. Many of these results were based on only a 
small fraction of spectrosc opic identification s of th e op- 
tically brightest o bjects. iGruppioni et al.l (|1999f ) and 
iCiliegi et~a l. (2003) have pointed out that, when identify- 
ing spectroscopically fainter objects, early-type galaxies 
constitute the majority. This is easily explained: if we 
cut our sample at V < 22, for example, we would be sam- 
pling a region of parameter space typical of SFG, with 
R ~ 0.5, z ~ 0.5, logP r - 23, and logL x - 41 (see 
Figs. [U [U and H]) . We would then wrongly infer that 
the sub-mJy population does not include the large num- 
ber of AGN, which become accessible only when probing 
fainter magnitudes. 

Our results are in broad ag reement with a num ber of 
recent papers. For example, ICiliegi et al.l |2003) have 
obtained a deep 5 GHz VLA image of the Lockman 
hole. The analysis of their 63 radio sources, which reach 
~ 50 fxJy, using a color-color diagram and a radio flux 
density - optical magnitude plot, suggests that AGN are 
st ill the dominant pop ulation at their flux density levels. 

ISmolcic et al.1 ()2008l ) have applied a new method to dis- 
entangle star forming from AGN galaxies using ~ 2,400 
1.4 GHz radio sources in the VLA-COSMOS survey. This 
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relies on a correlation between optical rest-frame colors 
of emission-line galaxies and their position in a spectro- 
scopic diagnostic diagram based on emission line flux ra- 
tios. They deduce that the radio population in the flux 
density range 50 [iJy —0.7 mJy is a mixture of 30 — 40% 
SFG, 50 - 60% AGN galaxies, and - 10% QSO, where 
the latter are defined as optical point sources. Their 
method works only up to z ~ 1.3 and therefore beyond 
t hat value their results are less secure. 

iSevmour et al.l (|2008f) have used a variety of (mostly) 
radio-based diagnostics, including radio morphology, 
spectral indices, radio to near-IR and mid-IR flux den- 
sity ratios, to discriminate on a statistical basis AGN 
from SFG in a deep VLA/MERLIN survey of the 13 ff 
XMM-Newton/ Chandra Deep Field, which includes 449 
sources. They find a dominance of SFG at the faintest 
flux densities (~ 50 ^Jy at 1.4 GHz), with AGN still 
making up around one quarter of the counts. They also 
find that the relative contribution of the two classes is 
flux density dependent, with SFG becoming the more 
numerous population below ~ 0.1 mJy. 

It is important to note that our selection criteria, by 
being not only radio-based, provide information on the 
overall nature of our sources, unli ke, for example, thos e 
of ISevmour et all (|2008l ) (see also iRichards et alj|2007h . 
which separate radio sources into those whose radio emis- 
sion is AGN and star formati on dominated. 

Regarding radio-quiet AGN. Uarvis fc Rawlingj (|2004D 
modeled the radio counts with a population of radio- 
quiet AGN becoming the dominant com ponent below 1 
mJy. We do not confirm this prediction. Simpson et al. 
(2006), on the other hand, suggested that > 20% of radio 
sources in the 0.1 — 0.3 mJy range were radio-quiet AGN 
and that these sources could be more important than 
SFG for the up-turn in the counts. While our results 
are consistent with the first point, they are not with the 
second one. 

We notice that radio-quiet AGN counts appear to be 
similar in shape to those of SFG, scaled down by a fac- 
tor ps 3 — 4 (see Fig. [7]), and quite different from those 
of the whole AGN population, which includes basically 
only radio-loud sources above 1 mJy. This argues against 
radio-quiet AGN being simply the faint end of a contin- 
uum of radio powers and should provide relevant infor- 
mation for the solution of the puzzle of the origin of their 
radio emission. 

6. SUMMARY AND CONCLUSIONS 

We have used a deep radio sample, which includes 256 
objects down to a 1.4 GHz flux density of 42 /zJy se- 
lected in the Chandra Deep Field South area, to study 
the nature of sub-mJy sources. Our unique set of ancil- 
lary data, which includes reliable optical/near-IR iden- 
tifications, optical morphological classification, redshift 
information, and X-ray detections or upper limits for a 
large fraction of our sources, has allowed us to shed new 
light on this long standing astrophysical problem. By 
analyzing the ratio of radio to optical luminosity and the 
radio and X-ray powers of the sources with morphologi- 
cal and redshift information, we have selected candidate 
star-forming galaxies and AGN from a complete sample 
of 194 objects. As the first two parameters by themselves 
are not very good discriminants between star-forming 
galaxies and AGN, optical morphology and especially X- 



ray data turned out to be vital in establishing the nature 
of faint radio sources. We have also complemented our 
data analysis with model calculations based on brighter 
samples and counts in other bands. Our main results can 
be summarized as follows: 

1. The well-known flattening of the radio number 
counts below 1 mJy is mostly due to star-forming 
galaxies, which are missing above ~ 2 mJy but be- 
come the dominant population below « 0.1 mJy. 

2. AGN exhibit the opposite behavior, as their counts 
drop at lower flux densities, going from 100% of 
the total at ~ 10 mJy down to < 50% at the 
survey limit. This is driven by the fall of radio- 
loud sources, as radio-quiet objects, which make 
up ~ 20% of sub-mJy sources, display relatively 
flat counts. 

3. Radio-quiet AGN make up about half of all AGN. 
Their counts appear to be a scaled down version, by 
a factor ps 3 — 4, of those of star-forming galaxies, 
and are very different from those of the radio-loud 
AGN population. This should provide a clue to the 
origin of their radio emission. 

4. Star-forming galaxies make up < 60% of sub-mJy 
sources down to the flux limit of this survey. This 
has to be regarded as a robust upper limit, as when- 
ever we had to make some assumptions, we chose 
to maximize their numbers. This result is at vari- 
ance with the many papers, which over the years 
have suggested a much larger dominance of these 
sources. On the other hand, our results are in 
broad agreement with a number of recent papers, 
which found a significant AGN component down to 
ps 50 /iJy. 

5. The results of our model calculations agree quite 
well with the observed number counts and pro- 
vide supporting evidence for the scenario described 
above. Moreover, they imply that sub-mJy radio- 
loud AGN are dominated by low-power, Fanaroff- 
Riley type I radio galaxies, as their high-power 
counterparts and radio-loud quasars are expected 
to disappear below ^0.5 — 1 mJy. 

6. We find a correlation between X-ray and radio 
power for star-forming galaxies, in agreement with 
most previous studies. 

7. Despite some previous claims, we find that the ra- 
tio of radio to optical luminosity depends more on 
radio luminosity, rather than being due to optical 
absorption, and is therefore not related to a high 
degree of extinction. 

Considering the apparent emerging population of low 
luminosity AGN at microjansky levels, care is needed 
when interpreting radio source counts in terms of the 
evolution of the star formation rate in the Universe. 

We plan to expand on this work by using our deeper 
(7 /xJy per beam ov er the whole Exten ded CDFS region) 
radio observations ((Miller et al. 200 8|) and the recently 
released 2 Msec Chandra data (|Luo et al.ll2008f h 
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